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We propose that the pseudogap state observed in the transition metal oxides can be explained 
by a three-dimensional flux state, which exhibits spontaneously generated currents in its ground 
state due to electron-electron correlations. We compare the energy of the flux state to other classes 
of mean field states, and find that it is stabilized over a wide range of t and 8. The signature of 
the state will be peaks in the neutron diffraction spectra, the location and intensity of which are 
presented. The dependence of the pseudogap in the optical conductivity is calculated based on the 
parameters in the model. 
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The motivation for this work is the observation of a 
pseudogap that opens up in optical conductivity mea- 
surements of the three-dimensional transition metal oxide 
SrRuOa [1] above its ferromagnetic transition tempera- 
ture of Tc ~ 150 K. A pseudogap has also recently been 
seen in BaRuOs [2]. In this pseudogapped regime, p (T) 
increases linearly with temperature, passing through the 
Ioffe-Regel limit without saturation [3], behavior indica- 
tive of a "bad metal" [4] . The optical conductivity in this 
state is proportional to the non-Fermi liquid behavior of 
uj^ 1 / 2 at high frequency and has a peak at low frequen- 
cies [1] at approximately 250 cm -1 , the precise location 
of the peak being temperature dependent. 

We propose that this pseudogap state can be under- 
stood by considering a ground state with spontaneously 
generated electronic currents circulating around the pla- 
quettes. The currents arise from electron-electron cor- 
relations, due to the bi-quadratic terms in the Hamil- 
tonian. The state which we propose is a generalization 
of the two-dimensional flux states invented by Affleck 
and Marston [5] , and studied in their chiral extension by 
Wen, Wilzcek, and Zee [6]. Unlike the two-dimensional 
case, there is no possibility of fractional statistics in 
three dimensions. However, the spontaneous generation 
of gauge fields is a possibility in three dimensions, and 
these gauge fields can lead to a ground state with cir- 
culating electronic currents. Earlier work was done on 
three-dimensional flux states by Laughlin and cowork- 
ers [7, 8] and Zee [9]. 

In actuality, SrRuOa has 5 bands crossing the Fermi 
surface formed by hybridizing the Ruthenium d orbitals 
with the Oxygen p orbitals [10]. The crystal structure is 
orthorhombic, becoming cubic at temperatures greater 
than 900 K [11]. Undoubtedly, the actual electronic 
structure of SrRuC>3, and particularly the presence of 
a van Hove singularity near the Fermi surface, influence 
the material's behavior. The model which we consider is 
vastly simplified and serves as a starting point for consid- 



ering the nature of the pseudogapped state in the three- 
dimensional transition metal oxides. A model which in- 
corporates some of these electronic features, but does not 
focus on the pseudogap regime, has been proposed by 
Laad and Muller-Hartmann [12]. 



I. MODEL SYSTEM 

The Hamiltonian that we consider is the single orbital 
t-J model given by 



a ]ai 



(1) 



where the sum over (ij) denotes nearest neighbors on a 
cubic three-dimensional lattice. Implicit in this equation 
is that we have set U — oo. The hopping matrix clement 
t which appears in Equation 1 is taken to be an effective 
hopping clement, which has been greatly reduced due to 
this on-site Coulomb repulsion. The value of t will be 
set by two calculations in this paper: the stability of 
the flux phase versus other mean field states calculated 
in Section II, and the value of the optical conductivity 
calculated in Section V. 

The spin operators may be written in terms of the 
fermion operators to give the Hamiltonian 

H = ~i E E c I c j*4« ,c ™' -'EE 4 c *r + ^n- 



(ij) aa' 



(ij) a 



The Hamiltonian will be treated in the mean field, or 
Hartree-Fock approximation. We make the replacement: 



E C ia c ja 



Xij 
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a J a 



Xij 



(3) 
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The assumption is made that the term in brackets, which 
corresponds to fluctuations about the mean-field Xij> is 
small and can be included only to linear order. The re- 
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FIG. 1: Unit cell for symmetry-broken bond Hamiltonian. 
The dark line passes through the four atoms in the planar 
unit cell. The spheres at the center of the cubes represent 
oppositely charged Dirac monopoles with the "tail" running 
through the interface between the two cubes. 



suiting Hamiltonian is given by 



H = \ E {l^l 2 - 2 [(*« + j) c h + H - c -] } 



(4) 



where we have dropped the redundant spin index. There 
is no a priori reason to believe that the fluctuations about 
the mean field will be small, although it has been rigor- 
ously shown in the two-dimensional case for the large n 
limit [5], where n is the particle spin. 

We allow the Xij to break the translational symme- 
try of the lattice. We choose a four-atom unit cell as 
shown in Figure 1. The lattice is generated by the prim- 
itive translation vectors (1,0,1), (1,0,-1), and (0,2,0), 
in units of the bondlength b. The Xij m Equation 4 arc 
then parameterized by 12 complex numbers. We use the 
notation Xw-, where the index i = 1, . . . , 4 is the location 
of the atom in the unit cell and the index v gives the 
direction {x, y, z}. 

This choice is made because it allows for the forma- 
tion of a 7r-per plaquette flux phase, something that a 
two atom unit cell does not allow in three dimensions. 
The model is computationally simpler than the eight- 
band model studied by Laughlin [7, 8], at the expense of 
picking out a preferred direction. To get a feel for what 
this corresponds to, one can think of the gauge fields in 
the sample as being generated by Dirac monopoles of al- 
ternating charge sitting at the centers of each cube. The 
"tails" of the monopoles are connected to form dipoles. 
The if direction in our model corresponds to the dipolar 
axis. 

The self-consistency of the model is the requirement 
that the energy as determined by Equation 4 be a local 
minimum with respect to variations of the twelve com- 
plex parameters Xiv This can be seen by writing 



self-consistency relation 



c \ C i+f 



(6) 



Note also that at this minimum, the expectation value of 
the Hamiltonian is given by 



Kf j 4 3 

w = -tEE 



| Xw | 2 + 2-^Re( Xw ) 



(7) 



The Hamiltonian in Equation 4 contains both a field 
strength and a term which we define as the bond Hamil- 
tonian: 



Hb = J E E (** + j) C Uja + H-C. 
(ij) ° V 



(8) 



The bond Hamiltonian can be diagonalized by introduc- 
ing a set of operators 



(9) 



where Af is the number of sites in the lattice, A is the 
band index which runs from 1 to 4 and the Uq\ (£) are a 
set of functions periodic in the unit cell. The bands are 
determined by the eigenvalue equation 



HqUqX = e qA WqA, 



(10) 



where £f q is written in terms of the hopping elements 
and Uq\ is a 4 component vector Uq\ (1) • • • Uq\ (4). For 
the unit cell depicted in Figure 1 , the Hamiltonian takes 
the form 



Hq = J 



( o m % o 

Vi o V4 
v* 3 o V2 



(ii) 



where 
Vi 

V2 



X 2y e 



-iq y 



X3 V e lqv + X* 4y e lq y 



Vi = X2xe^+X2ze iq '+xl x e- iqx +xl z e-^.(12) 

The tilde in the expression above means that these num- 
bers include the actual hopping clement: Xiv = Xiv+tj J ■ 
The eigenenergies for the four bands are given by 



(H) = 



4 6 

EE{ 



Xi, 



c i c l+ u ) +H.c. 



(5) 



Minimizing this function with respect to x* v gi yes the 



±J 



i 



± 



\vim-vam\ 



(13) 



with the usual Einstein summation convention. In Sec- 
tion III we will consider the band structure and eigen- 
states of one particular solution of the bond Hamiltonian, 
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the flux state where the hopping elements are complex 
and equal in magnitude: \x-iu\ = X f° r ai l h v - 

For t = the gauge fields \ij are unobservable. How- 
ever, as soon as t is increased they lead to real circulating 
electronic currents. Consider the site i and the six sites 
j which are its nearest neighbors. The sum of the cur- 
rents flowing outward from site i is then the time rate of 
change of the number operator m on site i 



0.25 



\dt 



2i 



7 -E( J ^'+ i )( c ^)- H - c -' ( 14 ) 



where the factor of 2 arises from the spin. The term 
Xij ( c ] c i) i s rea l an d gi yes no contribution, but for t > 
and Xij complex, the bonds carry electronic currents. 
These currents are the signature of the flux state. They 
have long-ange order and can be probed by neutron 
diffraction as will be shown in Section IV. 



II. THE PHASE DIAGRAM 



0.15 




0.05 - 



FIG. 2: Phase Diagram for Bond States. The diagram shows 
the competition between the four types of states considered 
in the text. The contours in the flux phase region show the 
average flux per plaquette, as defined by Equation 16. 



The Hamiltonian described in Section I admits a num- 
ber of self-consistent solutions. Wc have performed a 
numerical search as a function of the hopping element t 
and the doping S. The function to be minimized is 



E = — 



Xiv 



kA 



(15) 



where N m = W/4 is the number of unit cells in the lat- 
tice. The search is performed using Powell's method in 
the space of the \%v ■ At each point in the space, the bands 
are determined by Equation 13, and the lower magnetic 
bands are filled up with N e — Af (1 — S) electrons. 

Our search is limited to three classes of states. We 
consider the flux phase with \xi\ — \x\ f° r ai l i with 
the phases of the hopping elements unconstrained. At 
t = S = the flux phase has flux <& = tt per plaquette as 
defined by 



n 



Xi 



plaquette 



\Xi 



(16) 



Note that a flux of tt and a flux of — 7r per plaquette are 
indistinguishable since in cither case the electron acquires 
a phase of —1 upon traversing the plaquette. This is the 
chiral symmetry of the model. Away from t = S = 0, the 
flux per plaquette is no longer tt and the chiral symmetry 
is broken. 

Another state that we have considered is the dimer 
state. In this case each site forms a bond with a neigh- 
boring site. In the case of t = 5 = 0, one particular 
manifestation of the state is Xiz — Xix — 1 with all other 
Xiv = 0. The dimer state has flat bands e q = ± J and can 



be considered to be a charge-density state with the charge 
localized on the bonds for which Xiv = 1 [5] . Away from 
the point t = 5 = 0, the state becomes partially dimcr- 
ized. There are in principle a number of non-equivalent 
dimer solutions, but the solution mentioned above ap- 
pears to be the lowest energy configuration. 

At values of 6 > 0.2 another local minimum is the kite 
state which has xiy = X2x = X3z = Xiy ^ Xiv- The 
state is called a kite state because the lines of charge, 
if the analogy from the dimer state above is used, form 
zigzagging patterns through the lattice, and could lead 
to a lattice distortion from the large Coulomb repulsion 
between charged lines. This particular instantiation of 
the kite state is chosen for the same reasons as the dimer 
above. The last state which is considered is the uniform 
state with all Xiv equal and real valued. This state is 
a simple Fermi liquid, with a renormalizcd value of the 
hopping matrix element. 

The phase diagram is shown in Figure 2. While we 
have performed the calculation for an arbitrary flux state 
as described above, the phase diagram shows a calcula- 
tion performed using a restricted parameter set. The 
parameters used are shown as different arrows in Fig- 
ure 3 below. This is done because at S w 0.1, there is a 
transition to a flux phase with a different type of order- 
ing which is outside the scope of the current discussion. 
Apart from this transition, there is no qualitative change 
in the phase diagram when the unconstrained flux state 
is considered. From the phase diagram we see that the 
flux state is stabilized over a fairly wide range of dop- 
ing, but that the flux per plaquette decreases from the 
value of tt as one leaves the point t = 6 = 0. In the 
calculations which follow we will assume that t/J = 0.1, 
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a point at which the highly symmetric 7r-per plaquette 
flux state discussed below in Section III is a reasonable 
approximation to the actual mean field state. 

One must also consider the possibilities of other types 
of order which are not described by the mean-field 
The most insidious of these is antifcrromagnctic order. 
At t = 5 = 0, the energies of the two stabilized states 
discussed above are 



F --^1 

-^dimer ^ ; 



AfJ 

£ flux «-0.95— . (17) 



For comparison, the Neel state, which is characterized by 

= (-l)\ (18) 



Xij 







o c. 



has an energy of — 3A/"J/4, lower than either of the two 
bond states at t = 6 = 0. In order for the bond states to 
be actualized, a term has to be added which will desta- 
bilize the antiferromagnetic order. This can be done by 
adding a next-nearest neighbor hopping term J' [13], in 
which case the energy of the Neel state will be 



3AfJ 



1 - 2 



(19) 



The energies of the bond states are actually unchanged 
up to a value of J'/ J w 1/3, which is the threshold for 
acquiring a nonzero value of the next nearest neighbor 
Xij as shown by Laughlin and Zou [14]. At the value 
J' = J/3, the energy of the Neel state is equal to the 
energy of the dimer state and very close to the energy 
of the flux state. Doping will also serve to destabilize 
the antiferromagnet so that the crossover will actually 
occur at a lower value of J' than the one reported here. 
Therefore, while it is not treated explicitly in this paper, 
some term such as the next-nearest neighbor J' must be 
added to this model in order to make the bond states 
energetically favorable with respect to the Neel ordered 
state. 



III. THE FLUX PHASE 

In the calculations which follow, we consider the tt per 
plaquette flux phase described in Section II which is only 
truly stabilized at 6 = t = 0. Some care must be used 
in selecting this state since at t = 6 = the system is 
invariant under a gauge transformation where 



Xij 



(20) 



Away from this point, this symmetry disappears, and the 
low-energy state is the one for which the quantity 



J2 Re 



(21) 



is a maximum as can be seen from Equation 7. The state 
is shown in Figure 3. It can be found by either maximiz- 
ing the function in Equation 21 for an arbitrary gauge 




FIG. 3: The -k per plaquette flux phase. All bonds have the 
same magnitude. The arrows correspond to complex phases of 
4> (>) , 4> — it/A (>), and 37r/4 — <f> (►), where 4> is determined 
by tan 4> = \/2. This diagram shows the same portion of the 
lattice as is shown in Figure 1. 



transformation, or by numerically continuing a state from 
t > down to t = 0. The choice of the correct symmetry- 
breaking gauge is important since it will affect the dis- 
tribution of electronic currents in the sample and hence 
observable features such as the neutron diffraction and 
optical conductivity. 

In this state, the Hamiltonian in Equation 9 can be 
rewritten using a set of Dirac matrices a x , 
such that 



a yi and a z 



Hq = 2J Ixl cos ( qub "> au > 



where the matrices satisfy the algebra 



(22) 



(23) 



Explicit forms for the matrices are given in Appendix. It 
is found numerically that |x| ~ 0.4. The eigenvalues of 
the Hamiltonian are twofold degenerate and are given by 



e q = ±2J| X 




COS' 



{q v b). 



(24) 



The band structure is shown in Figure 4. At half-filling 
the Fermi surface reduces to two isolated points at qb = 
(7r/2,7r/2,7r/2) and qb = (tt/2, tt/2, -tt/2) shown at the 
point T, in Figure 4. The low energy excitations about 
these points are relativistic. 

In order to calculate the neutron cross-section and the 
optical conductivity, we also will need the eigenvectors 
which appear in Equation 9. The matrix u whose rows 
correspond to the bands 1 ... 4 and whose columns cor- 
respond to the position in the unit cell is given by 



u a = 



1 

71 



/Kl 


p-i4> 





i^i 7q \ 


l7ql 





33_ e ~i<t> 

l7ql 


-n J^3- 


M 


g— i<t> 







V l7ql 





_ 7q p -i<t> 





(25) 
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FIG. 4: Flux Phase band structure. Points in the Brillouin 
zone are T = (0,0,0), A = (tt/2, 0,0), R = (tt/2, 0, tt/2), 
and £ = (7r/2, 7r/2, 7r/2). There are two distinct Dirac points 
located at qb = (tt/2, tt/2, tt/2) and qb = (tt/2, tt/2, -tt/2). 



where we have defined the quantities 



»7q 



cos q y b 



7q =e 



. cosq x b — icos g z 6 



r q = 



cos 2 (g„6). 



(26) 



The phases of the eigenvectors in Equation 25 have been 
selected so that the eigenstates in Equation 9 are invari- 
ant under q — > q + Q where Q is any vector in the recip- 
rocal lattice. 



IV. NEUTRON SCATTERING 

The flux states can be probed by neutron scattering, 
as the neutron spin interacts with the magnetic dipoles 
generated by the real electron currents circulating on the 
plaquettes. The interaction potential [15] is written 



V (r) = 2 1 (4 c j - H - c -) ex P T i 3 A ' dl 



m 

where the vector potential is given by 



/X X 



|3 ' 



eh 
7 S. 



In these expressions S is the neutron spin and 7 
is a constant. It can be shown [17] that 



(27) 

(28) 
1.91 



/ 



dre iq ' r »y(r) 



(7ro) © 



(29) 



with the current operator defined as 

3 " = E c k+ q c k.^ ( k > q) 

k 

/„ (k, q) = cos (fc„6) - cos {k v b + q v b) . 



(30) 



In these expressions we have replaced the hopping ele- 
ment t by ti 2 /2m*b 2 , with m* the effective mass. Since 
we are assuming (see Section II) that t/J = 0.1, the effec- 
tive mass will be quite large. If we assume that J takes 
the typical value of 0.1 eV, we have m* w 15 m. The 
vector v is a unit vector with the sum running over the 
x, y, and z directions. 

Converting this to a cross-section and averaging over 
the spin, assuming that the incoming beam is unpolar- 
ized, the expression for the cross-section is given by 



da 
dtt 



(7ro) 2 (^V 
\m* / 



|q&l 



4 / j 



q v ' 



qv_ 



(31) 

In order to calculate the matrix elements in Equation 31, 
we need to rewrite the current operator in terms of the 
eigenstates of our system. We first break up the mo- 
mentum sum so that it runs only over the first Brillouin 
zone 



j » = E E c k+q+Q, c k+ Ql u (k + Qi, q) , 



(32) 



k i=l 



where Q16 = 0, Q 2 6 = (ir, 0, ir), Q36 = (0, n, 0), and 
Q46 = 7r. We can rewrite the electron operators at mo- 
mentum k + Q in terms of the c q (£) defined as 



4 



iq.(R+ r/ ) t 



C R+r t ' 



(33) 



R 



where R runs over all the unit cells and vt is the position 
of the £ th atom in the unit cell. This introduces a matrix 
g with elements gn — exp [iQi ■ re] /2, and results in the 
current operator 

5irra5ip c k+ q ( m ) c k (p) fv (k + Qi, q) • (34) 

k imp 

We also note that f v (k + Q j , q) = [Q v ] ■■ f v (k, q) , where 
the matrix Q v is diagonal with elements exp [iQi ■ ub}. 
The sum on i can then be performed to obtain 

J " = E E C U q H mp c k (p) /„ (k, q) . (35) 

k mp 

The matrix gQ^g in the above expression has a natural 
interpretation. It merely connects all sites in the lattice 
which are connected by a hopping element in the v di- 
rection. It can be written in the form 



gQ v g 



Finally, we can rewrite Equation 35 by inverting the 
eigenvector matrix in Equation 25. This gives 



J - = EEh {gQ v g) 4 +q 



k w 



yx .fAK q)^A,k+q^V,k- 

(37) 
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If we assume that the lower bands are completely filled, 
we can then write a simple expression for the expectation 
values of these matrix elements. 



W = L L Tr ' Uk ^Q v g) «l+q (k. q) V (38) 



The prime on the trace indicates that it runs only over 
the lower two bands, and the sum on r runs over all 
vectors in the reciprocal lattice. 

Expression 38 is quite general and can be used to cal- 
culate properties away from zero doping by restricting 
the sum on k such that k < kp. Additionally, this equa- 
tion assumes nothing about the actual structure of the 
eigenvector matrix Uk- Specializing to the case of the 
flux state discussed in Section III, the trace for v = y is 



1 



Tk+c 
7 k 



ry k sin ( 



(39) 



Here we have used the fact that rfc+q = — ?7k, a condition 
enforced by the presence of f y in Equation 38. We see 
that the trace is zero unless 7k+q = — 7k > which means 
that ( J y ) vanishes unless q v b is an odd multiple of tt for 
all v. For the v = x and v = z cases, the traces are the 
same as can be seen from Equation 36. They are given 

by 



?7k+q 



Vk J 



) Ini7k 



COS ( 



> — i ( 1 + k+q ) Re7k sin < 
Vk ) 



where we have taken 7k+ q = — 7k for the same reasons 
as above. This expression takes two different values de- 
pending on whether q y b is an even or an odd multiple of 

7T. 

Scattering will only occur at the reciprocal lattice vec- 
tors, as guaranteed by the delta function in Equation 38. 
The reciprocal lattice is shown in Figure 5. There is no 
magnetic scattering at the nuclear locations, since they 
occur at even multiples of ir. There is additionally no 
scattering at the line centers qb = try, since all three 
matrix elements vanish at these points. Scattering does 
occur at both the face-centers at q& = (tt, 0, tt) and at 
the body-center at qb = tt. The former distinguishes 
the scattering from Bragg scattering from a cubic anti- 
ferromagnet. It is likely that the actual material would 
consist of domains containing all dipolar orientations, so 
that scattering would actually be observed at all the face- 
centers. 

The magnitude of the scattering will in general be quite 
small. We write the cross section as 



da 



v 0m \m* / ^ 



(40) 



(41) 

where we have rewritten the Kronecker delta function 
from Equation 38 in terms of the Dirac delta function 
with the proper normalization of (27r) 3 /V. The term 
vom = 46 3 is the volume of the unit cell. The structure 
factor is given by 



M(q) = 

3 | q& |< 



(3 + cos q y b) 



1 



1 



(q x b) (q z b) z 



r 



1 — cos q y b 
(Qybf 



(1 — cos q x b) (1 — cos q z b) 



(42) 



In deriving this expression, we have repeatedly used the 
fact that 



1 ^ cos k v cos k v i 



k 



r k 



(43) 



with the Tk defined as in Equation 26. This relation 
follows from the symmetry of the momentum sum and 
Equation 15. The structure factor in Equation 42 takes 
the same value at the two smallest scattering angles cor- 
responding to the points Q 2 6 = (tt, 0, tt) and Q46 = tv in 
the reciprocal lattice: 



M(q) 



64 
3^4 



I XI 



(44) 



To give a feel for the order of magnitude of the scattering 
from the flux states, we compare to the scattering from 



an antifcrromagnct. In that case one has 



da 
dtt 



AF 



AF 



3 



V AF 



( 7 r ) 2 (^) 2 ^|f (r)\ 2 S(q-T) 

(45) 

where the factor of 2/3 arises by assuming that the sub- 
lattice magnetization is along a crystallographic axis. We 
can estimate the form factor for the antiferromagnet by 
assuming it is the same as that of chromium, a typi- 
cal band antifcrromagnct. Chromium has a form factor 
of F(iv/b) w 0.4 [16]. The unit cell in the flux phase 
is twice as large as the antiferromagnetic unit cell and 
therefore N^ F = 2Af m and Vq F = vo m /2. Assuming 
that the spins are 50% polarized such that (S v ) = 1/4, 
and taking m* = 15m, we see that the scattering from 
the antiferromagnet is roughly 170 times larger than the 
scattering from the flux state at the wave- vector qb = tt. 
Note that the discrepancy in size is due primarily to the 
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FIG. 5: Reciprocal Lattice Vectors for Flux Phase. The black 
dots show the scattering from the nuclear centers. The white 
dot shows the antiferromagnetic scattering vector at qb — 
7T. The shaded dots at the face centers show those points in 
the reciprocal lattice space of the flux phase which produce 
scattering, and the lighter shaded dots along the line centers 
show those points in the reciprocal lattice space which do not 
produce scattering. 



size of the effective hopping element t, which must be 
small compared to J if the flux phase is going to be sta- 
bilized away from 5 = as was shown in Section II. 



V. OPTICAL CONDUCTIVITY 

In our model, the peak in the optical conductivity 
arises from transitions between the bands shown in Fig- 
ure 4. The model is too simple to accurately predict the 
optical conductivity of a material such as SrRufJ3. How- 
ever, the calculation illustrates both the dependence of 
the location of the peak on the spin-exchange energy J 
and the intensity of the peak on the ratio of t/J. 

The calculation of the optical conductivity is very sim- 
ilar to that of the neutron diffraction. In this case we 
couple the system to a time-dependent vector potential 
A = A (i) i>, where the vector potential's time depen- 
dence and relation to the electric field are given by 



A(t) = Ae 



-iujt 



E 



ILO 



A. 



(46) 



We assume that the wavelength of the light is long enough 
that we can ignore any spatial dependence in the fields. 
The vector potential couples to the hopping terms in the 
Hamiltonian. A phase is acquired according to 



CjCj exp 



l 4- I A-d£ 

he 



(47) 



It is important to note that the corresponding Xij a P~ 
pearing in Equation 4 also acquire an equal and opposite 
phase so that when we expand the Hamiltonian to lin- 
ear order in the vector potential A (t) we do not get any 
contribution from the terms proportional to Xij- The 



perturbation to the Hamiltonian is 
H' 



(48) 



where L = A / ' 1/ ' 3 6 is the length of the sample and the 
current operator in the u direction is given by 



2iteb ^-^ 
3v = HlX^ 



[4 c i+v 



(49) 



Note that the operator in Equation 49 defines the total 
current, not the current density, flowing in the v direc- 
tion. The complex optical conductivity is related to the 
induced current density in the sample by 



Jind = 0"E. 



(50) 



The induced current can be calculated using linear re- 
sponse theory. In that case it can be shown from Equa- 
tions 46, 48 and 50 that the real portion of the opti- 
cal conductivity is related to the complex portion of the 
current-current correlation function 



CTi (w) = — Imx(w) , 



(51) 



where the current-current correlation function is given by 

1 



n 



hw - E n + E + is 



1 



Two + E n — E Q + is 



(52) 



The infinitesimal s arises from assuming that the pertur- 
bation in Equation 48 vanishes at t = — oo. The states 
|0) and \n) are to be evaluated at t = — oo or A = 0. The 
current operator in Equation 49 can be written in terms 
of the eigenstates of the system, following the same ap- 
proach applied in Equations 32 through 38 . The result 
is that 



J, = 4 



teb 
Lh 



[ ?ik (dQ V 9) u 

k AA' 



A' A 



sin {k v b) V'lkV'A'k- 



(53) 

Equation 53 is very similar to Equation 37. In this case, 
however, we are considering matrix elements connecting 
the ground state to excited states and so pick up the con- 
tributions at A ^ A'. Averaging the optical conductivity 
over the three directions in the lattice, it takes the form 



(7l (to) 



12\x\ 2hb \Jj M-M 



(54) 

where T q was defined in Equation 26 and where we 
have defined a dimensionless frequency given by \i — 
fou)/4:J\x\- The structure factor is given by: 

2 4 , 2 
v A=1A'=3 
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sition metal oxides. Its signature will be the presence of 
weak neutron diffraction peaks arising from the ordered 
electronic currents in the material. The fact that SrRuC>3 
is not near an antiferromagnetic transition and the fact 
that the three dimensional flux states produce scattering 
at wavevectors other than qb = ir make the system an 
excellent candidate in which to observe this type of or- 
der. Further theoretical work is warranted to understand 
how the actual electronic structure of the material will 
influence the behaviors discussed here. 



FIG. 6: Interband contribution to optical conductivity for the 
flux phase. 
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The sums in Equation 55 can be performed to give 
^(q) - (^+ (7k ~ 6 ^ )2 ) (sin 2 fc^ + sin 2 fc z 6) 



+ 3 (! + »£) sin 2 k y b, 



(56) 



with the quantities defined as in Equation 26. 

Equations 54 and 56 have been evaluated numerically 
and the result is shown in Figure 6. Within our model, 
the location of the peak is proportional to the exchange 
energy J. Assuming this takes the value 0.1 eV, we find 
that the peak occurs at approximately 1000 cm -1 , which 
gives order of magnitude agreement with the observed 
value of 250 cm -1 in SrRuC>3 [1]. Due to the simplic- 
ity of the model we are solving, one would not expect 
more accurate agreement. The magnitude of the peak 
is governed by the ratio of t/J. Direct comparison of 
this quantity with experiment is more difficult. This is 
due to the fact that SrRu03 has five bands, whereas we 
have considered only a single orbital. Additionally, our 
calculation only considers the interband contribution to 
the conductivity, whereas the real material also has an 
intraband contribution from thermally excited carriers. 
If one assumes that the result shown in Figure 6 needs 
to be scaled by a factor of roughly five to account for 
the number of orbitals in SrRuC>3, the results are reason- 
able compared with the measured conductivity of 6000 
JJ-'cm- 1 . 
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APPENDIX 

The matrices from Equations 22, 23, are given explic- 
itly by 



( Ze~ 3iw / 4 \ 

Z*e™/ 4 

Z * e 3iw/4 

V Ze-™/ 4 J 

( Ze-™/ 4 \ 

z*e 3 ™/ 4 

z * e m/4 



a x 



a z = 



a v = 



\ Ze- 31 "/ 4 

/ Z \ 

Z* 

Z* ' 

\ Z / 











(A.l) 



where 



VI. DISCUSSION 

We propose that the three-dimensional flux state is a 
good candidate for the pseudogap state seen in the tran- 



Z = e 40 = -4= 



1 

V3 



(A.2) 
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